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Abstract 

In this article, we recalculate the contributions from the vacuum condensates up to dimension- 
£C) t 6, including the one-loop corrections to the quark condensates a s {qq) and partial one-loop 

corrections to the four-quark condensates a s (qq) , in the operator product expansion, then 
study the masses and decay constants of the pseudoscalar mesons D, D a , B and B s with the 
CN \ QCD sum rules. The predictions f D = (208± 11) MeV and f B = (189± 15) MeV are consistent 

with the experimental data within uncertainties, while the prediction fo s = (241 ± 12) MeV 
is still below the lower bound of the experimental data fo a ~ (260.0 ± 5.4) MeV. 
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Q> 1 Introduction 

^ ', The charged pseudoscalar mesons ir + , K + , Df and B + mesons can decay to a charged lepton pair 

r"| ■ t + vg through a virtual W + boson. To the lowest order, the decay width is 
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V{P^lv) = ^f pm jm P [l-^ \V qiq2 \\ (1) 

where the mp is the P meson mass, the f P is the corresponding decay constant, the is the 
' £ mass, the V qiq2 is the Cabibbo-Kobayashi-Maskawa matrix element between the constituent 

quarks 91(721 an d the Gp is the Fermi coupling constant. The CLEO collaboration obtains the 
values f D = (202 ± 41 ± 17) MeV [J, (222.6 ± 16.7±1|) MeV [2], (205.8 ± 8.5 ± 2.5) MeV [5J from 
the decay D+ -> fj+^; f Ds = (259.5 ± 6.6 ± 3.1) MeV, JdJId = 1.26 ± 0.06 ± 0.02 [4] from the 
decay D+ -> fi+v^; f Ds = (252.5 ± 11.1 ± 5.2) MeV [5], (259.7 ± 7.8 ± 3.4) MeV [6J from the decay 
£>+ -> t+v t . The Babar collaboration obtains the value f Ds = (258.6 ± 6.4 ± 7.5) MeV [7 from 
the decays D~ — > i~vi. The Belle collaboration obtains the value fp> s = (275 ± 16 ± 12 MeV [5] 
from the decay Df Now the average values listed in the Review of Particle Physics are 

f D = (206.7 ±8.9) MeV, f Ds = (260.0 ± 5.4) MeV and f D Jf D = 1.26 ±0.06 [9]. 

There have been many theoretical works on the decay constants of the heavy pseodoscalar 
mesons, such as the QCD sum rules (QCDSR) [TUl [III [TH UHl HH H3 HBl H3 HB] , the lattice QCD 
(LQCD) [H HB [23 the Bethe-Salpeter equation (BSE) [23 H, the relativistic potential 
model (RPM) rjHHlllTj, the field-correlator method (FCM) [25], the light-front quark model 
(LFQM) [2S1 [30], the chiral extrapolation [3T], the extended chiral-quark model [32], e t c - There 
are discrepancies between the theoretical values (from QCDSR and LQCD) and the experimental 
data, which maybe signal some new physics beyond the standard model [15) . In the QCD sum 
rules, the expressions in the operator product expansion are different from each other in one way 
or the other, this maybe due to different approximations are taken by the authors [TT1 [UH [T51 [55] . 
In this article, we recalculate the contributions from the vacuum condensates up to dimension-6, 
including the one-loop corrections to the quark condensates a s {qq) and partial one-loop corrections 
to the four-quark condensates a 2 s (qq) 2 , in the operator product expansion, then study the masses 
and decay constants of the pseudoscalar mesons D, D s , B and B s with the QCD sum rules. The 
QCD sum rules is a powerful theoretical tool in studying the ground state hadrons [541, [55] , 

The article is arranged as follows: we derive the QCD sum rules for the masses and decay 
constants of the heavy pseudoscalar mesons in Sect. 2; in Sect. 3, we present the numerical results 
and discussions; and Sect. 4 is reserved for our conclusions. 
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2 QCD sum rules for the heavy pseudoscalar mesons 



In the following, we write down the two-point correlation functions II(p) in the QCD sum rules, 

H(p) = iJd 4 xe l P x {0\T{j 5 (x)4(0)}\0), (2) 

Js(x) = Q(x)i"f5q(x) , (3) 

where the pseudoscalar currents Js(x) interpolate the heavy pseudoscalar mesons, Q = c,b and 
q = u,d, s. We can insert a complete set of intermediate hadronic states with the same quantum 
numbers as the current operators J§{x) into the correlation functions Tl(p) to obtain the hadronic 
representation [M| [33] . After isolating the ground state contributions from the heavy pseudoscalar 
mesons, we get the following result, 

n (p) = ? 5V + "-. ( 4 ) 

y ' (m Q + m q ) 2 {m 2 p -p 2 ) K ' 

where the decay constants fp are defined by 

(o\Mo)\p(p)) = fp f p ■ (5) 

Now, we briefly outline the operator product expansion for the correlation functions Tl(p) 
in perturbative QCD, and use the charmed-strange (or bottom-strange) mesons to illustrate the 
procedure. We contract the quark fields in the correlation functions Tl(p) with Wick theorem 
firstly, 

H(p) = i J d 4 xe^Tr{ To S l3 {x) To S%{-x)) , (6) 

where the Sij(x) and S^Ax) are the full quark propagators, and can be written as 

iSjj _ 5 l3 m s _ Sjj , iSjj j6m s {ss) _ 5jjX 2 (sg s aGs) iSjjX 2 #m s (sg s aGs) 
blAX > ~ 2tt 2 z 4 ~ 4ir 2 x 2 ~ 12 {SS) + 48 ~~ 192 + 1152 

32ir 2 x 2 3456 ' U 



<% 9sG n a ^ a»0Qi + m Q ) + Q6 + m Q )a^ 



* jW (2tt) 4 J \H-rriQ 4 {k 2 - m 2 Q ) 2 

g s D a G n 0x t^(f^ a + f Xa ?) g 2 a ^)ijG a ap GlM aP,lv + /"^ + Z"^) 1 
+ 3(fc 2 - m 2 Q Y A{k 2 - m 2 Q f + " ' J ' 

^ = (^m Q ) 7 A (^mQ)7^ + fflg )/(^mQ), 

= (# + m Q ) 7 ^ + mQ) 7 ^ + m Q ) 7 ^ + m Q ) 7 ^ + m Q ), (8) 

and t n = 4-, the A™ is the Gell-Mann matrix, the i, j are color indexes, D a = d a — ig s G™t n [33] ; 
then complete the integrals both in the coordinate space and in the momentum space; finally obtain 
the correlation functions II(p) at the level of the quark-gluon degrees of freedom. The analytical 
perturbative 0(a s ) corrections [10] and semi- analytical perturbative 0{a 2 ) corrections [12] to the 
perturbative term are available now. In this article, we take the elegant expression in Ref.[10]. For 
the vacuum condensates, we calculate the Feynman diagrams shown explicitly in Figs. 1-5 in the 
operator product expansion, where the solid and dashed lines represent the light and heavy quark 
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Figure 1: The perturbative 0(a s ) corrections to the quark condensate (ss). 



propagators, respectively. In calculation, we take the light quark mass m q as a small quantity 
and expand it perturbatively. In Fig.l, there exist divergences, the quark condensate in the full 
propagators should be replaced as 



<££> 
12 



12 



1+1 r 



(9) 



In this article, we take the dimension D = 4 — 2e to regular the divergences. In calculations, 
we observe that the mixed condensates (see Fig. 2) are depressed by additional powers of 1/T 2 
compared to the quark condensates (see Fig.l). The perturbative 0(a s ) corrections to the mixed 
condensates are doubly depressed by the factor a s /T 2 and play a minor important role, and they 
are neglected in this article. In the massless limit, the second Feynman diagram in Fig. 3 does 
not contribute to the gluon condensate ( "'^ ), the QCD spectral density is f ( ° sGG ) instead of 
^( a ^ G ). In calculating the fifth Feynman diagram in Fig. 3, we use the light quark propagator 
S(k) in the momentum space, 



S(k) 



(g 3 s GGG) 
48 



if6 



(10) 



The expression presented in Ref.[3S] is in four-dimension, we add the factor |e to obtain the 
propagator in /^-dimension as there are divergences. In calculating the Feynman diagrams in 
Figs. 4-5, we use the equation of motion, D v G a iLV = J2 q =u d s 9sQ1^ a Q} an d take the approximation 
(ss) = (qq), furthermore, we take assumptions of the vacuum saturation and factorization [34 , 
and use the following formulae, 



(qi^qq^q) = 



(11) 



The factor cannot be neglected when companied with divergences in the loop integral. In Fig. 5, 
we present the Feynman Diagram cannot be written as the perturbative 0(a s ) corrections to the 
four-quark condensates a s (ss) 2 (shown in Fig. 4). The four-quark condensates a s (ss) 2 play a minor 
important role, the perturbative 0{a s ) corrections a 2 (ss) 2 can be safely neglected, although they 
appear in one- loop order. 

Once analytical expressions of the correlation functions at the quark level are obtained, we 
can use the dispersion relation to obtain the QCD spectral densities, then take the quark-hadron 
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Figure 2: The diagrams contribute to the mixed condensate (sg s aGs). 








Figure 3: The diagrams contribute to the gluon condensate ( " sGG ) and three-gluon condensate 
{glGGG). 
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Figure 4: The diagrams contribute to the four-quark condensate (ss) 2 of the order 0(a s ). 




Figure 5: The typical diagram contributes to the four-quark condensate (ss) 2 of the order 0(a 2 s ). 



duality and perform the Borel transforms with respect to the variable P 2 
QCD sum rules, 
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(12) 



where 



2Li 2 (a;) + logxlog(l - a;) 



-log 

2 x 



log(l 



a; log 



-logic , 



(13) 



r(0,x) = e' x / dt e'\ 

Jo t + x 

f x 1 

Li 2 (x) = - / dt- log(l - i) , 



(14) 



and the sq is the continuum threshold parameter. The perturbative 0(a s ) corrections R(x) are 
taken from Ref.|10). We can derive Eq.(12) with respect to 1/T 2 , then eliminate the decay constant 
fu B to obtain the QCD sum rules for the mass m^ s . The QCD sum rules for the decay constants 
and masses of the pseudoscalar mesons D, B and B s can be obtained with simple replacements. 



3 Numerical results and discussions 

The masses of the pseudoscalar mesons listed in the Review of Particle Physics are m B ± = (1869. 5± 
0.4)MeV, m D o = (1864.91 ±0.17) MeV, m D ± = (1969.0 ± 1.4) MeV, m B ± = (5279.25 ± 0.26) MeV, 
m B o = (5279.55 ± 0.26) MeV, m B o = (5366.7 ± 0.4) MeV 0. In 2010, the Babar collaboration 
observed four excited charmed mesons L>(2550), L>(2600), £>(2750) and D(2760) in the decay 
modes J D°(2550) -> D*+ir~, £>°(2600) -)• D*+ir- , D+tt-, D°(2750) -» D*+n-, L>°(2760) -» 
D+tt-, D+(2600) -)■ D°tt+ and L>+(2760) £>°7r+ respectively in the inclusive e+e~ cc 
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interactions at the SLAC PEP-II asymmetric-energy collider [35]. The doublet (D(2550), D(2600)) 
are tentatively identified as the 2S doublet (0~, 1~) [37]. So we can take the threshold parameters 
as S£, = 6.2 GeV 2 and s^, = 7.3 GeV 2 tentatively to avoid the contaminations from the high 
resonances, here we have taken into account the width of the Z?(2550) and the SU(3) symmetry 

mo = 



breaking effects. If additional uncertainties 5sq = 0.5 GeV are supposed, then s 



(0.5 - 0.7) GeV and 




tod. 



(0.6 — 0.8) GeV, the contributions from the ground states are 



fully included. In Ref. 18 a , S. Narison takes the threshold parameters as s° D = (5.3 — 9.5) GeV 
and s B = (33 — 45) GeV 2 . In this article, we take threshold parameters as s° B — (33.5 ± 1.0) GeV 2 
and s° B ^ = (35.0 ± 1.0) GeV 2 for the bottom mesons, then - m B = (0.4 - 0.6) GeV and 

■> B — m Ba = (0.5 — 0.6) GeV. The contributions from the ground states are fully included and 
the contaminations from the high resonances are very small if there are some contaminations. We 
expect that the couplings of the pseudoscalar currents to the excited states are more weak than 
that to the ground state mesons. For example, the decay constants of the pseudoscalar mesons 
7r(140) and 7r(1800) have the hierarchy: / w (i3oo) -C /7r(i4o) from the Dyson-Schwinger equation 
[55] . the lattice QCD [35], the QCD sum rules [3D], etc, or from the experimental data [4"T] . 

The vacuum condensates are taken to be the standard values (qq) = — (0.25±0.01 GeV) 3 , (ss) = 



(0.8 ±0.1) (q>g), (qg s aGq) = ml{qq), (sg s aGs) = 



ss 



ml 



(0.8 ±0.1) GeV at the energy scale 



H = 1 GeV |42j . The quark condensate evolves with the renormalization group equation, (qq)(^ 2 



{m){Q 2 ) 



,{Q) 



The value of the gluon condensate 



asGG ) has been updated from time to time, 
we use the recently updated value (S^G) = (0.022 ± 0.004) GeV 4 03], 



and changes greatly 

and take the three-gluon condensate as (g 3 s GGG) = (8.8±5. 5) GeV 2 (a s GG) = (0.616±0.385) GeV 



In this article, we take the WS masses m c (m 2 ) = (1.275 ± 0.025) GeV, TO 6 (m 2 ) = (4.18 ± 
0.03) GeV, m s (ii = 2 GeV) = (0.095 ± 0.005) GeV from the Particle Data G roup [9], and set 
m q = 0. Furthermore, we take into account the energy-scale dependence of the MS masses from 
the renormalization group equation, 



m s (p 2 ) 
m c (/i 2 ) 



m s (4GeV 2 ) 
m c (m 2 ) 



m b (ml) 
1 



a s {2GeV) 

12 



a s (m c ) 



a s (m b ) 

h logt 6 2 (log 2 1 - log< - 1) + 6 6 2 



b 2 t 



b 4 t 2 



(15) 



where t = los 



b 



33-2n/ 

127T • 



153-19n/ 



2857- 



and 339 MeV for the flavors n/ 

/i = y m? D — m? c w 1 GeV; for the 6-quark, we take rif = 4 and /i = -J m B — m 

In Fig. 6, we plot the contributions from different terms in the operator product expansion with 
variations of the Borel parameters. From the figure, we can see that the convergence of the operator 
product expansion cannot be satisfied for the D and D s (B and B s ) mesons at the region T 2 < 
0.9 GeV 2 (T 2 < 3.0 GeV 2 ). In Figs. 7-8, we plot the masses and decay constants with variations of 
the Borel parameters at large ranges. Although there appear minimum platforms for the masses 
and decay constants of the D and D s mesons at T 2 < 0.9 GeV 2 , the Borel windows cannot be 
chosen in such regions. For the B and B s mesons, the decay constants decrease monotonously 
with increase of the Borel parameter at T 2 < 4 GeV 2 . We choose the suitable Borel parameters to 
satisfy the two criteria (pole dominance and convergence of the operator product expansion) of the 



■21, ~ ■ - i >8^ A = 213 MeV, 296 MeV 

5, 4 and 3, respectively [9J. For the c-quark, we take rif = 3 and 

2.5 GeV. 
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T 2 (GcV' z ) 


s (GeV") 


pole 


mp(GeV) 


/p(MeV) 


D 


1.2 - 1.8 


6.2 ±0.5 


(66 - 93)% 


1.87 ±0.10 


208 ± 11 


D s 


1.2 - 1.7 


7.3 ±0.5 


(80 - 96)% 


1.97 ±0.10 


241 ± 12 


B 


4.7-6.1 


35.5 ± 1.0 


(46 - 70)% 


5.28 ±0.06 


189 ±15 


B s 


5.1 - 6.6 


35.0 ± 1.0 


(47-71)% 


5.37 ±0.06 


216 ± 16 



Table 1: The Borel parameters, continuum threshold parameters, pole contributions, masses and 
decay constants for the heavy pseudoscalar mesons. 





/r> (Me V) 


f Da (MeV) 


/s(MeV) 


/s s (McV) 


Id 3 1 Id 


/s s / Ib 


Expt [§] 


206.7 ±8.9 


260.0 ±5.4 


194 ±9 




1.26 ±0.06 




QCDSR 14 


177 ±21 


205 ± 22 


178 ±14 


200 ± 14 


1.16 ±0.16 


1.12 ± 0.11 


QCDSR pi] 


206.2 ± 7.3 


245.3 ± 15.7 


193.4 ± 12.3 


232.5 ± 18.6 


1.193 ±0.025 


1.203 ±0.020 


QCDSR 18 


204 ±6 


246 ±6 


207 ±8 


234 ±5 


1.21 ±0.04 


1.14 ±0.03 


LQCD [20, 


197 ±9 


244 ±8 






1.24 ±0.03 




LQCD [21 


213 ±4 


248.0 ±2.5 


191 ± 9 


228 ± 10 


1.164 ±0.018 


1.188 ±0.018 


LQCD [22 


218.9 ± 11.3 


260.1 ± 10.8 


196.9 ±8.9 


242.0 ±9.5 


1.188 ±0.025 


1.229 ±0.026 


BSE 23 


238 


241 


193 


195 


1.01 


1.01 


RPM 26 


234 


268 


189 


218 


1.15 


1.15 


FCM [28 


210 ± 10 


260 ± 10 


182 ±8 


216 ±8 


1.24 ±0.03 


1.19 ±0.02 


LFQM 30, 


205.8 ±8.9 


264.5 ± 17.5 


204 ±31 


270.0 ±42.8 


1.29 ±0.07 


1.32 ±0.08 


This work 


208 ±11 


241 ± 12 


189 ±15 


216 ± 16 


1.16 ±0.07 


1.14 ±0.11 



Table 2: The decay constants of the heavy pseudoscalar mesons from the experimental data and 
some theoretical calculations. 

QCD sum rules, and reproduce the experimental values of the masses. The threshold parameters, 
Borel parameters, pole contributions and the resulting decay constants are shown explicitly in 
Table 1. 

In Table 2, we compare the predictions of the present work with the experimental data and 
some (not all) theoretical calculations. The value Jb — 194 ± 9 MeV listed in the Review of 
Particle Physics [9] is the average of the lattice QCD calculations [22j [45] ■ The present values 
fn = (208 ± 11) MeV and Jb = (189 ± 15) MeV are consistent with the experimental data within 
uncertainties, while the value fr> s = (241 ± 12) MeV is below the lower bound of the experimental 
data fr> a — (260.0 ± 5.4) MeV In the present case, JdJId ~ /bJIb, the heavy quark 
symmetry works well. Most of the theoretical predictions (including the present prediction) of the 
ratio fD s /fD are below the experimental data. In Ref.|16j. A. Khodjamirian estimates the upper 
bound fo < 230 MeV and fo 3 < 270 MeV based on the QCD sum rules, the present predictions 
satisfy the constraints. The existence of a charged Higgs boson or any other charged object beyond 
the standard model would modify the decay rates, for example, the leptonic decay widths are 
modified in two-Higgs-doublet models [35] ■ The discrepancies between the theoretical calculations 
and the experimental data maybe smeared, if there are new physics beyond the standard model 
and the theoretical estimations of the decay constants are robust. 

4 Conclusion 

In this article, we recalculate the contributions from the vacuum condensates up to dimension-6, 
including the one- loop corrections to the quark condensates a s (qq) and partial one-loop corrections 
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Figure 6: The contributions from different terms in the operator product expansion, the A, B, 
C and D correspond to the D, D s , B and B s , respectively; while the 1, 2, 3, 4, 5, 6 and 7 
denote the contributions from the perturbative terms, quark condensate, mixed quark condensate, 
gluon condensate, four-quark condensate 0(a s ), three-gluon condensate and four-quark condensate 
O{o 2 s ) 1 respectively. 
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Figure 7: The masses of the heavy pseudoscalar mesons, the A, B, C and D correspond to the 
D, D s , B and B s , respectively. 
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Figure 8: The decay constants of the heavy pseudoscalar mesons, the A, B, C and D correspond 
to the D, D s , B and B s , respectively. 
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to the four-quark condensates a 2 s {qq) 2 , in the operator product expansion, and obtain the analytical 
QCD spectral densities, then study the masses and decay constants of the heavy pseudoscalar 
mesons with the QCD sum rules. The values f D = (208 ± 11) MeV and f B = (189 ± 15) MeV are 
consistent with the experimental data within uncertainties, while the value fr> s — (241 ± 12) MeV 
is still below the lower bound of the experimental data fr> s = (260.0 ± 5.4) MeV. If the present 
predictions are robust, new physics beyond the standard model are need to smear the discrepancy. 



Acknowledgements 

This work is supported by National Natural Science Foundation, Grant Number 11075053, and 
the Fundamental Research Funds for the Central Universities. 



References 

[1] G. Bonvicini, Phys. Rev. D70 (2004) 112004. 
[2] M. Artuso et al, Phys. Rev. Lett. 95 (2005) 251801. 
[3] B. I. Eisenstein et al, Phys. Rev. D78 (2008) 052003. 
[4] J. P. Alexander, Phys. Rev. D79 (2009) 052001. 
[5] P. U. E. Onyisi et al, Phys. Rev. D79 (2009) 052002. 
[6] P. Naik et al, Phys. Rev. D80 (2009) 112004. 
[7] P. del Amo Sanchez et al, Phys. Rev. D82 (2010) 091103. 
[8] L. Widhalm et al, Phys. Rev. Lett. 100 (2008) 241801. 
[9] J. Beringer et al, Phys. Rev. D86 (2012) 010001. 
[10] T. M. Aliev and V. L. Eletsky, Sov. J. Nucl. Phys. 38 (1983) 936. 

[11] C. A. Domingucz and N. Paver, Phys. Lett. B197 (1987) 423; S. Narison, Phys. Lett. B198 
(1987) 104; L. J. Reinders, Phys. Rev. D38 (1988) 947; M. Jamin and M. Munz, Z. Phys. C60 
(1993) 569; S. Narison, Phys. Lett. B520 (2001) 115; M. Jamin and B. O. Lange, Phys. Rev. 
D65 (2002) 056005; H. Y. Jin, J. Zhang and Z. F. Zhang, Phys. Rev. D81 (2010) 054021. 

[12] K. G. Chetyrkin and M. Steinhauser, Phys. Lett. B502 (2001) 104; K. G. Chetyrkin and M. 
Steinhauser, Eur. Phys. J. C21 (2001) 319. 

[13] A. A. Pcnin and M. Steinhauser, Phys. Rev. D65 (2002) 054006. 

[14] J. Bordes, J. Penarrocha and K. Schilcher, JHEP 0412 (2004) 064; J. Bordes, J. Penarrocha 
and K. Schilcher, JHEP 0511 (2005) 014. 

[15] S. Narison, Phys. Lett. B668 (2008) 308. 

[16] A. Khodjamirian, Phys. Rev. D79 (2009) 031503. 

[17] W. Lucha, D. Melikhov and S. Simula, Phys. Lett. B701 (2011) 82; W. Lucha, D. Melikhov 
and S. Simula, J. Phys. G38 (2011) 105002. 

[18] S. Narison, larXiv: 1209.20231 



11 



[19] A. X. El-Khadra et al, Phys. Rev. D58 (1998) 014506; D. Becirevic et al, Phys. Rev. D60 
(1999) 074501; S. Collins et al, Phys. Rev. D63 (2001) 034505; K. C. Bowler et al, Nucl. Phys. 
B619 (2001) 507; G. M. de Divitiis et al, Nucl. Phys. B672 (2003) 372; A. Gray et al, Phys. 
Rev. Lett. 95 (2005) 212001; C. Bernard et al, PoS LATTICE2008 (2008) 278; E. T. Neil et 
al, PoS LATTICE2011 (2011) 320; D. Becirevic et al, JHEP 1202 (2012) 042. 



[20 
[21 

[22 
[23 
[24 
[25 
[26 
[27 
[28 
[29 
[30 
[31 
[32 

[33 
[34[ 
[35 
[36 
[37 

[38 
[39 
[40 
[41 
[42 

[43 
[44 



B. Blossier et al, JHEP 0907 (2009) 043. 

C. T. H. Davies et al, Phys. Rev. D82 (2010) 114504; H. Na et al, Phys. Rev. D86 (2012) 
034506. 

A. Bazavov et al, Phys. Rev. D85 (2012) 114506. 

Z. G. Wang, W. M. Yang and S. L. Wan, Nucl. Phys. A744 (2004) 156. 

G. Cvetic, C. S. Kim, G. L. Wang and W. Namgung, Phys. Lett. B596 (2004) 84. 
P. Colangelo, G. Nardulli and M. Pictroni, Phys. Rev. D43 (1991) 3002. 

D. Ebcrt, R. N. Faustov and V. O. Galkin, Phys. Lett. B635 (2006) 93. 
M. Z. Yang, Eur. Phys. J. C72 (2012) 1880. 

A. M. Badalian, B. L. G. Bakker and Yu. A. Simonov, Phys. Rev. D75 (2007) 116001. 

H. M. Choi, Phys. Rev. D75 (2007) 073016. 

C. W. Hwang, Phys. Rev. D81 (2010) 114024. 

X. H. Guo and M. H. Weng, Eur. Phys. J. C50 (2007) 63. 

D. Ebert, T. Feldmann, R. Friedrich and H. Reinhardt, Nucl. Phys. B434 (1995) 619; S. Nam, 
Phys. Rev. D85 (2012) 034019. 



A. Hayashigaki and K. Terasaki, hep-ph/0411285 



M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, Nucl. Phys. B147 (1979) 385. 
L. J. Reinders, H. Rubinstein and S. Yazaki, Phys. Rept. 127 (1985) 1. 
P. del Amo Sanchez et al, Phys. Rev. D82 (2010) 111101. 

Z. G. Wang, Phys. Rev. D83 (2011) 014009; B. Chen, L. Yuan and A. Zhang, Phys. Rev. 
D83 (2011) 114025; A. M. Badalian and B. L. G. Bakker, Phys. Rev. D84 (2011) 034006; P. 
Colangelo, F. De Fazio, F. Giannuzzi and S. Nicotri, Phys. Rev. D86 (2012) 054024. 

A. Hocll, A. Krassnigg, C. D. Roberts and S. V. Wright, Int. J. Mod. Phys. A20 (2005) 1778. 

C. McNeile and C. Michael, Phys. Lett. B642 (2006) 244. 

K. Maltman and J. Kambor, Phys. Rev. D65 (2002) 074013. 

M. Dichl and G. Hiller, JHEP 06 (2001) 067. 

P. Colangelo and A. Khodjamirian, |hep-ph/0010175j ; B. L. Ioffe, Prog. Part. Nucl. Phys. 56 
(2006) 232. 

S. Narison, Camb. Monogr. Part. Phys. Nucl. Phys. Cosmol. 17 (2002) 1. 

S. Narison, Phys. Lett. B693 (2010) 559; S. Narison, Phys. Lett. B706 (2012) 412; S. Narison, 
Phys. Lett. B707 (2012) 259. 



12 



[45] E. Gamiz et al, Phys. Rev. D80 (2009) 014503. 

[46] A. G. Akeroyd and C. H. Chen, Phys. Rev. D75 (2007) 075004; A. G. Akcroyd and F. 
Mahmoudi, JHEP 0904 (2009) 121. 



13 



